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Precision spectroscopy of trapped HfF + will be used in a search for the permanent electric dipole moment of the electron (eEDM). 



While this dipole moment has yet to be observed, various extensions to the standard model of particle physics (such as supersym- 
metry) predict values that are close to the current limit. We present extensive survey spectroscopy of 19 bands covering nearly 
' ^ ' 5000 cm 1 using both frequency-comb and single-frequency laser velocity-modulation spectroscopy. We obtain high-precision 
Q^rovibrational constants for eight electronic states including those that will be necessary for state preparation and readout in an 
i actual eEDM experiment. 
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1. Introduction 
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' r/l The existence of a permanent electric dipole moment of an 
>^»electron (eEDM) would have profound implications for funda- 
i-£h mental physics since it violates parity and time reversal sym- 
i ^ 'merries 01 Q [Hi- The current experimental limit of about 
1 x 10~ 27 e x cm 0, 0] is many orders of magnitude larger than 
the predicted eEDM from the standard model fl, but exten- 
sions to the standard model (such as supersymmetry) predict a 
dipole moment on the order of 10~ 29 to 10~ 26 e x cm jg]. This 
means that measurement of the eEDM provides a rigorous test 
of extensions to the standard model without having to disen- 
tangle new results from standard model predictions. HfF + and 
C^> the related ThF + have been suggested as candidate species in 
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such a measurement using trapped molecular ions |7[ 
,— ' The low-lying, metastable 3 Aj states in HfF + and ThF + have 
high sensitivity to the eEDM due to the large effective electric 
J^j field felt by one of the unpaired electrons when the molecule is 
polarized in low laboratory electric fields. As initially shown 
C3 by DeMille et al for PbO, the level structure of Q. = 1 states 
provides valuable checks to reduce systematic errors 01 ll . The 
use of molecular ions makes trapping straightforward, thus en- 
abling measurements with long coherence times. 

The general approach for the eEDM experiment with HfF + 
or ThF + will be to perform a Ramsey-type spectroscopy mea- 
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surement of the energy separation of two \mf\ — 3/2 magnetic 
sublevels of the 3 Ai 7=1 level (see Figure [TJ- In applied 
magnetic and electric fields these sublevels are split both by the 
Zeeman effect and by the interaction of the dipole moment of 
the electron (assumed to be parallel to the electron spin) with 
the internal electric field of the ion. This magnitude of this shift 
is given by 2d L ,& mt (using the convention of |9j] instead of that in 
jlOtl), where d e is the eEDM to be measured and £,■„, is the ef- 
fective internal electric field experienced by the electron (about 
24 GeV/cm Jgtl). HfF + will initially be formed and trapped in 
the X'Z + state; then optical (and possibly microwave) transi- 
tions will be used to create a coherent superpostion in the two 
desired nif sublevels of the 3 Ai state. After a free-evolution 
period a readout pulse will be applied to meaure the accumu- 
lated phase difference between the two sublevels, which is pro- 
portional to the energy difference between the sublevels. A de- 
tailed analysis of the proposed experiment can be found in IToTl . 
For this scheme, optical transitions from the metastable 3 Ai 
state will most likely be required for state-selective readout via 
laser induced fluorescence or resonant multiphoton photodisso- 
ciation (states (b) or (c) in Figure [TJ respectively); in addition, 
transitions coupling the X'E + state with the 3 Aj state via some 
excited state (state (a) in Figure [TJ will be necessary for state 
preparation in the 3 Aj state. 

Prior to this work, very little spectroscopic information was 
available for either HfF + or ThF + . The energies of the low- 
lying 3 Ai, 3 A2, and 3 A3 states were recently measured using 
pulsed-field ionization zero kinetic energy (PFI-ZEKE) spec- 
troscopy 1 12 , T3I1. Additionally, several isoelectronic species 
including TiO IllfHfl^l, TiF + and TiCP EEll have been 
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Figure 1 : Energy levels of HfF + relevant to the eEDM experiment. The 

eEDM experiment will utilize the J = 1 level of the 3 Ai metastable state. Mixed- 
spin excited states that couple well to both the ground 'Z + and the 3 Ai could 
be used for population transfer (state (a)). Readout from the 3 Ai could be ac- 
complished directly via laser-induced fluorescence (b) or resonantly-enhanced 
multiphoton photodissociation (c). The inset shows the structure of the / = 1, 
F = 3/2 manifold of the 3 Ai level in applied electric and magnetic fields. The 
eEDM experiment will measure the difference in energy between the upper 
(and lower) pairs of \mf \ = 3/2 states. The applied magnetic field splits these 
states by 3g/iB, where B is the magnitude of the magnetic field, n is the Bohr 
magneton, and g is the magnetic g-factor. In addition, the electron experiences 
an effective internal electric field (£,,«), which results in an additional splitting 
with magnitude 2t/ f ,£,„ f , where d e is the eEDM (aligned with the electron spin). 



studied. Importantly though, there are no fully characterized 
optical transitions in either HfF + or ThF + at wavelengths below 
1 //m (above 10 4 cm -1 ) where laser-induced fluorescence detec- 
tion is possible. Due to the large number of electrons involved, 
quantum calculations are challenging: current high-level calcu- 
lations have errors of perhaps 1000 cm -1 Ir9l [l9ll . 

The combined lack of experimental data and large theo- 
retical uncertainties necessitated precision spectroscopy over 
a very broad spectral range. To achieve this, we recently 
developed frequency-comb velocity-modulation spectroscopy 
(comb-vms), which provides simultaneously broad bandwidth, 
high resolution, high sensitivity, and ion specificity lEoll . We 
were able to measure 1000 crrT 1 of continuous spectra around 
12000 crrT 1 using comb-vms and observed four bands of 
HfF + . This information, combined with PFI-ZEKE data and 
theoretical predictions, was then used to guide specifically 
targeted scans based on single-frequency (cw-laser) velocity- 
modulation spectroscopy, which resulted in 15 additional bands 
with origins ranging from 9950 to 14600 crrT 1 . We obtained 
precise molecular constants for the 3 Aj, 3 Ilo-, 3 Ilo+, 'IT, 
3 Ili, 3 £ 0+ , and 3 <E>2 states, many of which will be of use for the 



eEDM experiment. We were also able to measure the electronic 
contributions to isotope splittings. In addition, we discuss the 
observed A-doublings, including the first measurement of the 
size of the small doubling in the 3 Ai. This value is important 
for determining the necessary applied electric field to fully po- 
larize HfF + as well as for estimating other systematic effects in 
the final eEDM experiment. Finally, the large amount of data 
enabled refinements in the ab initio calculations; we briefly dis- 
cuss the modifications and compare the new theory with the 
measurements. This may help to improve the accuracy of other 
calculations involving high-Z atoms. 

Outside of the eEDM experiment, precision spectroscopy of 
molecular ions is useful in a wide range of fields including fun- 
damental physics, chemistry, and astrophysics. In astrochem- 
istry for example, at least 22 ions (both positive and negative) 
ranging in complexity from Hi to H2COI-F and have been 



identified in interstellar and circumstellar gases [21]. Many of 
these ions are highly reactive and are believed to be interme- 
diates in a variety of reactions; nonetheless, some rate con- 
stants and branching ratios are still not well known ll22ll . For 
example, the pathway for formation of a very simple molecu- 
lar ion, CH + , remains elusive 12311 . Searches for new species, 
including efforts to identify the origin of the diffuse interstel- 
lar bands [24, 25b, benefit from laboratory measurements of 
optical transitions as well as measurements of rotational con- 
On the 



28] 



stants to aid microwave spectroscopy II26L 127 , 
physical chemistry side, precision spectroscopy of above 
the barrier to linearity (near 9913 cm -1 ) provides rigorous tests 
for ab initio theory [29, [3(1 31]. Many other carbocations 



are interesting for both their presence as intermediates in reac- 
tions such as combustion as well as the challenges they present 
to theory; for example, the spectrum of CH^ , a highly non- 
classical carbocation, still remains unassigned due to both its 
complexity and spectral interference from other contaminant 



species 1321 133L 13411 . A combination of comb-vms with a re- 
cently developed ion-beam spectrometer capable of sensitive, 
sub-Doppler vms of rotationally cold species ( 113511 ) could en- 
able mass-selective spectroscopy of ions over a broad band- 
width with many potential applications. 

2. Experiments 

2.1. Comb-vms 

Comb-vms combines cavity-enhanced direct frequency 
comb spectro scopy 11371,1381, 13911 with velocity-modulation spec- 



4311 for discrimination between ions and 



troscopy S HTf 
neutral species. This technique could also be readily extended 
to the detection of radicals by concentration modulation. By re- 
solving every comb mode of a femtosecond Tksapphire laser 
(3 GHz mode spacing) simultaneously over a wide spectral 
bandwidth, we have high resolution and absolute frequency ac- 
curacy with over 1500 channels measured at once. Figure |2] 
provides an overview of the experimental setup, more details 
can be found in lEoll . Briefly, we couple light from the comb 
into a ring cavity containing an ac discharge cell, which forms 
the HfF + and modulates the ions' drift velocity. The resulting 
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Figure 2: Setup for comb-vms. We use a Tiisapphire comb with a repetition 
rate of 3 GHz, which is transferred via fiber to the optical cavity. The comb 
is coupled into one direction of a ring bowtie cavity, which contains an ac- 
discharge cell with Brewster-angled end windows. The cavity finesse is about 
100 and the length is matched to be an integer multiple (25) of the comb cav- 
ity. We use liquid crystal variable retarders and a polarizing beamsplitter to 
alternate the direction of propagation through the cavity. The ac-discharge is 
driven at 10 kHz with about 500 mA. We flow helium gas through the tube at a 
pressure of about 6 Ton. HfF + ions are created by heating about 0.5 g of HfF4 
powder to 550 °C inside the discharge tube. The cavity transmitted light, which 
contains the modulated absorption signal, is analyzed using a two-dimensional 
imaging system and a lock-in camera. A cw Ti:sapphire laser is used to stabi- 
lize the comb and serves as a frequency reference in the spectra. More details 
can be found in (2oTl . A liber-coupled grating spectrometer is used to monitor 
the discharge fluorescence. 



modulated Doppler shift produces an absorption signal that is 
modulated at the discharge frequency. The comb light transmit- 
ted through the cavity is then spectrally dispersed using a high- 
resolution (1 GHz), two-dimensional cross-dispersion system 
I 44 ], which resolves every comb mode. We then image this 
onto a lock-in camera fl45[ |46fl to demodulate the absorption 
signal simultaneously on many comb teeth. 

We produce HfF + by heating about 0.5 g of HfF4 powder 
to about 550 °C in a 1 m long home-built discharge cell. He- 
lium buffer gas is flowed through the discharge such that the 
total pressure is about 6 torr. By recording the emission of the 
discharge using a low-resolution grating spectrometer, we were 
able to reliably monitor molecule production. With each load- 
ing of the oven we were able to run for about 3 hours. A 2.5 m 
long bow-tie optical enhancement cavity consisting of two 
flat, 98% mirrors (input and output couplers) and two 100 cm 
radius-of-curvature, low-dispersion, 99.9% high reflectors sur- 
rounds the discharge cell. The reflectivity of the input and out- 
put couplers was chosen to match the losses from the Brewster- 
angled windows and thus provide efficient input coupling. The 
cavity length is actively stabilized to an integer multiple (25) 
of the frequency-comb laser cavity length, which ensures that 
each comb component is coupled to a cavity mode. Our feed- 




Figure 3: Setup for cw-vms. We use a cw Ti:sapphire laser that is tunable 
from about 9500 - 14800 cirT 1 (1050 - 675 nm). After fiber coupling between 
tables, the laser is split with a 50/50 non-polarising beamsplitter. Each beam 
goes through an additional 50/50 non-polarizing beamsplitter (to separate the 
return light) and is then sent in opposite directions through the discharge cell. 
Half of the returning light in each direction is reflected by the beamsplitters 
and is differenced using an auto-balanced photodetector 13611 . This difference 
signal is then sent to lock-in detectors to record signal in-phase and 90° out-of- 
phase with the HfF + absorption. We placed irises in the beam paths to reduce 
the amount of discharge fluorescence observed on the auto-balanced detector; 
this helped both to reduce noise and primarily to reduce baseline drifts in the 
spectra. 



back loop provides about 20 kHz of bandwidth by using a high- 
speed, low-range piezo-electric transducer (PzT) ll47ll and a sec- 
ond long-range PzT. We use liquid-crystal variable retarders 
and a polarizing beamsplittler to rapidly (50 ms) switch the di- 
rection of propagation through the discharge tube while main- 
taining the comb-cavity lock. This allows us to subtract out 
slowly varying noise due to drifts in camera pixel offsets and 
also to improve the rejection of neutral background absorption. 

We record about 150 cirT 1 of spectrum spread over 1500 
channels simultaneously. One comb tooth is locked to a sta- 
ble cw Ti: sapphire laser, which provides an absolute frequency 
reference in our spectrum. For each measurement, we average 
and subtract images for each direction of propagation and also 
record the power per comb tooth by applying a calibrated am- 
plitude modulation to the laser. In order to fully sample the 
spectrum, we interleave 30 measurements with the cw laser 
stepped over 3 GHz. This results in a spectrum that covers 
150 crrT 1 sampled every 100 MHz with an absolute frequency 
accuracy of 30 MHz (set by a rubidium referenced waveme- 
ter). When the wavemeter was not calibrated for a particular 
measurement, we set the uncertainty at 100 MHz. For strong 
bands with fully resolved isotope structure, absolute accuracy 
of the determination of the band origin is limited predomi- 
nantly not by statistical errors but by the absolute knowledge of 
the cw Ti:sapphire laser frequency we use as a reference. For 
our purposes we were satisfied with the 30-100 MHz (0.001 - 
0.003 cirT 1 ) absolute accuracy of the wavemeter. This could 
be readily improved if needed. For analysis purposes we then 
interpolate the spectrum onto a fixed 0.001 cirT 1 grid. This al- 
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Figure 4: Comparison of coherent subtraction and single direction mea- 
surements. The signal to noise for coherent subtraction (using a 52/48 beam 
splitter) surpasses that of single direction measurements by a factor of 10 when 
the noise is primarily technical light intensity noise. Both measurements of a 
single Nj line were made with the cw laser. The dashed lines are a calculation 
of the expected lineshape (see text). The modification in the lineshape for the 
coherent subtraction is due to a stray differential phase shift of approximately 
7t/20 between the two counter-propagating beams, which results in a sensitivity 
to both absorption and dispersion. 

lows us to easily average or combine a collection of different 
scans. One full scan takes about 30 minutes and results in a 
single-pass fractional absorption sensitivity of 3 x 10~ 7 . Since 
one scan contains 45000 channels, this equates to a sensitivity 
of 4 x 10~ 8 Hz~'/ 2 (spectral element) -1 ^ 2 , which is the sensitiv- 
ity that a single-frequency laser system would need, in addition 
to being able to scan 150 crrT 1 continuously, to match the per- 
formance of the comb-vms system. We recorded spectra over 
1000 cirT 1 with both the oven on and off to check for contami- 
nation. 

2.2. Single-frequency vms 

The information obtained from the comb scans, combined 
with theory predictions, allowed us to scan other bands using 
single-frequency velocity-modulation spectroscopy as sketched 
in Figure [3] For these measurements we removed the en- 
hancement cavity and counter-propogated beams from the cw 
Ti: sapphire laser through the discharge tube. We then sub- 
tract these two beams using an auto-balanced photodetector 
l36Tl . Due to noise from the discharge, we reached a sensitiv- 
ity of about 5 x 10~ 8 Hz 1/2 with 1 mW on the detector; this 
noise seemed to be related to optical pickup from the discharge 
emission, which was reduced using irises, and to non-common- 
mode amplitude noise from acoustical pickup. Using this tech- 
nique we were able to find 15 more bands ranging from 9950 to 
14600 cirT 1 without having to scan the full spectral range. We 
can continuously scan about 0.5 crrT 1 in 4 minutes, which is at 
least 30 times slower than the comb-vms system. 

2.3. Coherent Subtraction 

We have also investigated a novel technique for differential 
detection that relies on coherent interference between the two 



counter-propogating laser beams instead of subtraction of pho- 
tocurrents as discussed previously. This is accomplished by 
coherently splitting and recombining the two beams using one 
non-polarizing beamsplitter to form a Sagnac-type interferome- 
ter containing the discharge cell and cavity if desired (blue inset 
to Figure 3]i. The signal is then detected at the destructive in- 
terference, or "dark port", of the interferometer. In the regime 
where laser intensity noise dominates, the signal-to-noise using 
coherent subtraction increases as the splitting ratio approaches 
50:50 between the two beams until detector readout noise or 
shot noise dominates. Physically, the fractional signal is in- 
creasing while the fractional noise level remains constant and 
dc power are both decreasing. We have tested coherent subtrac- 
tion versus a single direction of propagation using a photodi- 
ode and single-channel lock-in detector with our cw laser and 
demonstrate a factor of 10 gain in the signal to noise (blue data 
in Figure|4]i, which is the predicted improvement with the 52/48 
beamsplitter that was used. 

Coherent subtraction, unlike subtraction of photocurrents, is 
sensitive to phase as well as amplitude of the light field; thus, 
the resulting lineshape is dependent on both absorption (ap- 
proximately first-derivative shaped) and dispersion (approxi- 
mately second-derivative) and can vary if there are differential 
phase shifts between the two directions. We performed a sim- 
ple simulation, shown as dashed lines in Fig. [4] that reproduces 
the lineshape modification. These simulations were done by ap- 
plying a cosine modulation in time to the center frequency for 
both the absorption (assumed Gaussian lineshape) and disper- 
sion terms (obtained via Kramers-Kronig) and then selecting 
the appropriate frequency term from the Fourier cosine trans- 
formation. Repeating this at a variety of simulated laser fre- 
quencies generates the lineshape. The amplitude and linewidth 
of the single-direction simulation was scaled to the measured 
value and these values were used for the coherent-subtraction 
simulation. The constant phase offset (most likely due to the 
beamsplitter) added to the coherent subtraction simulation was 
chosen to fit the measured lineshape (about n/20 radians). In 
addition to the phase sensitivity, which could be useful for some 
applications, coherent subtraction has several other advantages. 
First, the dark port can be used to reduce the optical power seen 
by the detector thus avoiding saturation while still maintaining 
a shot-noise limited signal-to-noise ratio. Second, the subtrac- 
tion is performed prior to the detector and is wavelength in- 
dependent. This means that coherent subtraction can be used 
with comb-vms or other dispersive detection systems for in- 
creased signal-to-noise. We did not implement coherent sub- 
traction with the comb-vms system here because we were close 
to camera-noise limited, thus the gains would not have been 
very significant. 

3. Results and Discussion 

Figure |5ja) shows the spectrum of HfF + acquired using the 
comb-vms system in blue (offset for clarity) as well as all data 
obtained with single-frequency vms in red. As illustrated in 
the zoomed region shown in the lower panel of Figure |5jb), 
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Figure 5: All data obtained with cw- and comb-vms (red and blue, respectively). All y-axes are fractional single-pass absorption (xlO~ 5 ) and all x-axes are in 
cm . In (a) we show all of the HfF + spectra acquired to date. The comb-vms data is vertically offset for clarity. The lower panel of (b) shows a zoomed region 
of the comb data where some lines were also measured with the cw system. The upper panel of (b) plots predicted line positions from the fits; different colors 
correspond to different bands (magenta, «— 'S + (0,1); green, 3 rii «— 'Z + (3, 1); cyan, 3 n i <— 'S + (2, 0)) while different symbols correspond to Hf isotopes: 
-o, I80 Hf; -X, I79 Hf; -*, I78 Hf; -+, I77 Hf; — , 176 Hf. Inset (c) shows the 3 <D 2 «- 3 A] (0, 0) band (on the left) and part of the 3 X Q+ <- 'S+ (1, 0) on the right. The 
region in (d) illustrates the cancellation of isotope shifts in Av = transitions; the sharpness of P- and Q-branch lines is due to a cancellation between the rotational, 
vibrational, and electronic contributions to the isotope shift. The splitting in the (higher-/") R-branch lines that is not observed in the (lower-/") P- and Q-branch 
lines is due to A-doubling in the 3 Ai state. Inset (e) shows a prominent band-head from the 3 <I>2 <— 3 Ai (0, 0) band as well as weaker lines from the 3 £~ + «— 'S + (1, 
0) band. 
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the spectrum from about 12100-12300 crrT 1 is extremely con- 
gested due to the presence of many bands, each with five iso- 
topes, and the high temperature of our oven, which results in 
observed /" values up to about 70. The dynamic range of 
frequency-comb velocity-modulation spectroscopy is demon- 
strated by our ability to identify the overlapping 1 1 <— 'Z + 
(v' = 0, v" = 1) and 3 Ili <- [ E + (3,1) bands despite the dif- 
ference in linestrengths and an offset in band origin of only ~1 
cm -1 . Since 150 cirT 1 sections are acquired simultaneously 
when using the comb, relative linestrengths within the region 
are not influenced by variability in oven and discharge condi- 
tions, which significantly helps to disentangle bands. The up- 
per panel of Figure[5lb) shows predicted line positions from fits 
to three different bands (the fitting is discussed below), which 
illustrates our ability to resolve each isotope (— o, 180 Hf; -x, 
179 Hf; -*, 178 Hf; -+, 177 Hf; — , 176 Hf) for multiple bands. Fig- 
ure |5Jb) also demonstrates excellent overlap between the comb 
and cw spectra. 

Since we were fitting a variety of different transitions with 
various combinations of Hund's cases, we chose to fit each band 
with the general expression given in Eq. [TJ 

v(v' , v" , / , J" ,s',s") = v + F v (J',s')~ F v » (J",s") ( 1 ) 

where the rotational energies, F v are 

F V (J, s) = (B v - *|)/(J + 1) - (D v - s k -^)J 2 (J + If. (2) 

We include rotation (B), centrifugal distortion (D), A-doubling 
(k), and distortion of the A-doubling (ko) for each vibrational 
level, v; s is an e/f-symmetry term where s = +1 for e- 
symmetry and -1 for f-symmetry. k and ku are generic constants 
for the A-doubling since we have fit states with several differ- 
ent values of A and Q. In So states, only one symmetry term 
is chosen and k and ko are both set to 0. vo includes both the 
electronic and vibrational energies: 

vo = T' e - T'J + E vib , (3) 

E vib * K(v' + 1/2) - aj e x' e (v' + 1/2) 2 ) 

- K'(v" + 1/2) - cj e x' e '(v" + 1/2) 2 ). (4) 

Here T e is defined as the energy difference between the min- 
imum of the potential curve from the minimum of the X'S + 
curve and oj e and io e x e are the usual vibrational constant and 
anharmonic correction, respectively. 

Due to the congested nature of the spectra, fits were per- 
formed to a manually identified list of isolated lines (position 
only) for each band and for isotope in bands with well resolved 
isotope structure. Figure |6jb) shows the residuals from a rep- 
resentative fit; root-mean-square values for the residuals were 
typically < 0.005 cirT 1 with no apparent structure. For most 
Av = transitions the fits were performed to estimated isotope 
averaged line centers and then the rotational constants were 
shifted by reduced mass to the 180 Hf values. The 3 <I>2 <— 3 Ai 
bands exhibited line doubling at high-J due to A-doubling (dis- 
cussed in detail below), for these bands the center of the doublet 



was used for the fits. Similarly, the 1 II i < — 3 (0, 0) band exhib- 
ited doubling of every transition (i.e., 6 branches) due to large 
A-doubling in the 'IT; for this band all six branches were fit 
simultaneously by assuming the splitting to each transition was 
k' x ]'(]' + 1). A summary of the fitted constants for each band 
is given in Table Q] 

Assignments of the observed bands were not particularly 
straightforward, partially because the presence or absence of 
low-J lines usually could not be determined due to the high tem- 
perature and complicating isotope structure. In fact, only for the 
'ill <— 'S + (0, 0) band near 13000 cirT 1 were we able to directly 
establish that Q' = 1 and £2" = since we observed an R(0) 
line but no P(l) line. We also measured several bands with no 
Q-branch, which were assigned as 3 IIo+ <— 'X + and 3 E~ + <— 
This information, combined with data from ||l2ll . isotope shifts, 
and A-doubling structure provided the assignments in Table [TJ 
Due to possible ambiguities, several bands are still unassigned 
as indicated. 

From the fitted constants for the 16 assigned bands, we can 
determine constants for the X'S + , 3 Ai, 3 n -, 3 n +, 'IT, 3 H U 
3 Sq + , and 3 ®2 states as given in Table [2] 7q and AG1/2 are 
directly obtained from the measurements: Tq is defined as the 
energy of the v = level of an excited state relative to the X'Z + 
v = level, and AG1/2 is the energy difference between the 
v = and v = 1 levels. All other (equilibrium) constants are 
extracted from the data. The rotational constant for a given 
vibrational level is given to first order in v as B v — B e — a e (v + 
1 /2), where B e is the equilibrium rotational constant. Since 
we have measured at least two vibrational levels for each state, 
we can determine a e ~ B v - B v +i> an d then use this to obtain 
B e for each state. T e , a> e , and <jj e x e are determined using Eqs. 
[3] and [4] For the 3 IIi state, we were able to calculate u> e and 
u> e x e directly; for all other states we assumed a Morse potential 
so that (jj e x e - a 2 e <±> 2 e / '36B 3 + a e a> e /3B e + B e (i.e., the Pekeris 
relationship). 

In addition, we have observed isotope shifts of the state ori- 
gins, called 6T e in Table|2] that we attribute to Hf electronic iso- 
tope shifts (due to the finite nuclear charge radius) ll48[ l49i l50ll . 
We define the isotope shift between 180 HfF + and 178 HfF + for a 
given spectral line as 



= (1 - —)E ro , + (1 - J^)E vih + 5T e - ST'J. (5) 

Ml78 V ^ 178 

Here jUiso and jUns are the reduced masses, E ro , is the rotational 
energy in 180 HfF + (i.e., offset from the band origin), and E v n, 
is the vibrational energy in 180 HfF + . For bands with well re- 
solved isotope structure (i.e., Av + 0), ST' e - ST" was deter- 
mined either by fitting both vi 80 and v ( ' 78 or by measuring 5E 
near the origin and then subtracting the vibrational contribu- 
tion. For Av = bands, we found the frequency where the 
lines were the sharpest, indicating that 6E * (as illustrated 
by the Q- and P-branch lines in Figure Eld)); this location rel- 
ative to the band origin is used for E rot in Eq. [5] which then 
gave 5T' e - ST'J. Values of \5T' e - 6T' e '\ were < 0.1 cirT 1 for 
all transitions. To obtain the 6T e values given in Table |2j we 
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Figure 6: Residuals for three fits to the 3 IIo- <— 3 Ai (1, 0) band. The residuals are only plotted for assigned lines in each branch: R-branch (+), Q-branch(square), 
and P-branch (*). (a) Fit with no A-doubling terms included. Notice the large structure to the residuals, (b) Only k" included, i.e., only A-doubling proportional to 
J 2 in the lower state. No structure in the residuals; the scatter is primarily due to error in identification of the line center, (c) Only k'^ included, i.e., only A-doubling 
proportional to J 4 in the lower state. There is some slight structure to the residuals, indicating that the primary contribution to the A-doubling goes as J 2 . 



set ST"(X l I, + ) = and calculated relative shifts of the other 
states. The simplistic model of Eq. [5] neglects vibrational- 
band-specific perturbations or other more complicated level- 
dependent effects; however, results were consistent across dif- 
ferent transitions, which supports the assignment of the effect 
to electronic states, not individual bands. Except for the X'E + 
state the isotope shifts are all fairly similar with an average 
value of about -0.06 cm _1 (-1.7 GHz), which agrees well with 
the value for the Hf + 5d 2 6s 2 -> 5d 2 6s6p of -1.8 GHz Q. 

Due to our ability to see high-/ lines, we were able to mea- 
sure precise values for the A-doubling, parametrized by k and 
ko in Eq. [2] in several different electronic levels. In Hund's case 
(a) k can be related to the more standard A-doubling parameters 
defined in 115 111 for n states (o, p, and q) and to 5 a as defined in 
1 5211 for 3 A states; k was not measurable in <t> states. The val- 
ues of k and ko were useful in assigning many transitions and 
also provide some insight into inter-state interactions. For the 
eEDM measurement, the most important A-doubling parameter 
is the splitting of the / = 1 levels in the 3 Ai, as this goes into 
determining several different potential systematic effects. In a 
case (a) 3 Ai state, the A-doubling is expected to be given by 
±<5a x J (J + 1), where the upper and lower signs refer to the 
two e/f-symmetry levels 115 211 . We can make a rough estimate 
the value of 5 a by a simple scaling from that measured in TiO 
orWC JH: 



isotope splitting. In this band, the lines are not sufficiently well 
resolved to permit an accurate determination of the A-doubling 
constant, but by fitting the 3 no- <— 3 Ai (1,0) transition, we were 
able to precisely determine the sign and size of the A-doubling. 
As shown in Figure [6] the best fit was obtained with k" instead 
of k'p, implying that the A-doubling goes mainly as /(/ + 1) as 
expected (adding k'^ to the fit with k" resulted in no improve- 
ment). This fit gives a value of k" = -1.23(6) x 10~ 5 cm -1 , 
which gives o A = k" /2 = 6.2(3) x 10~ 6 cnr 1 = 185(9) kHz. 
The precision on 5a is only about an order of magnitude worse 
than that obtained for TiO using pure rotational spectroscopy 
1 1511 . illustrating the ability to obtain high-resolution spectro- 
scopic information from broad-bandwidth ro-vibronic spectra. 

A-doubling in the 3 no can arise from interactions with both 
'S + and 3 S + states, which result in a splitting of the O = 0~ and 
+ levels that is independent of J. We measured this splitting to 
be 189.099(4) cirT 1 with the + (e-symmetry) higher in energy 
than 0~ (f). To the extent that the X'S + state is a pure (scr) 2 
configuration, we can assume that this splitting is from a higher- 
lying 3 E + state with an {scrdcr) configuration. The spin-orbit 
interaction between the 3 Sq_ (f) and the 3 no- can be estimated 
by 



(i<<n (J5)n(J3)\\al + s-\a, (J3)cr 2 (a))) 2 l^ 5d) 



AE 



AE 



(7) 



OA 



£ 2 (Hf 2+ )fi 2 (HfF + ) 
<r 2 (WC/TiO)fi 2 (WC/TiO) 



x 5a(WC/TiO). (6) We have used (n(P)\l + s lo^a)) « V6 in the pure-precession 



This scaling relationship arises because the oa term results from 
the application of two spin-orbit operators (which scale as the 
atomic spin-orbit, and two L-uncoupling operators, which 



model 15 1L I55H . Using the measured splitting and the 
atomic spin-orbit coefficient, we estimate AE m 11000 cm -1 , 
which gives an approximate location of the 



3 V + 



Eq_ state 



to 



scale as the rotational constant B 1521 15311 : however, it neglects 
changes in state order and spacing. Using this, with £ obtained 
from |54], we predict 5 a ~ 50 kHz. The A-doubling in the 
3 Ai can be observed at high-J in the 3 <I>2 *~ 3 Ai (0, 0) band by 
comparing neighboring P-, Q-, and R-branch lines, as shown in 
Figure |5jd). Since the R-branch lines have the highest J, there 
is noticeable doubling in these lines compared to the P- or Q- 
branch lines, which indicates that this doubling is not due to 



be 21200 cm -1 . This is close to the predicted position of 
21694 cuT 1 from the new calculations (see Section |4j. In a 
Hund's case (a) basis, this splitting corresponds to the more fa- 
miliar 2(o + p + q) ifsTh. 

The cause of the A-doubling in the 3 ni state is primarily 
interaction with the 3 Sq_ state through the 3 no. This type of 



interaction is characterized by the parameter p in 15111 . With e- 
symmetry above f in 3 no and 3 no located below 3 ni, we would 
expect e above f in 3 ni, which corresponds to a negative value 
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of k as observed. The magnitude of k is expected to be about 
4Bp/An, where An is the spin-orbit parameter for the 3 n mani- 
fold and 



p * CCn^BL^lX^aL-S^Uoy/iEn - £ s ) (8) 



+ |3t 



CB({5d)(7T + \l + \cr 2 ) 2 /(En - 



(9) 



Here, C is a numerical factor dependent on the spin 05 IB . 
Again using the pure-precession model, we estimate k ~ 8 x 
10 -4 cm -1 , in reasonable agreement with experiment. 

We attribute most of the A-doubling in the to interac- 
tion with the nearly-degenerate nominal 3 X„~ + state. It is rea- 
sonable that the 1 El and 3 E~ are not pure (scrdn) and (d6) 2 
configurations, which would lead to increased spin-orbit in- 
teraction between 'ill and 3 2 Q+ . Estimating the magnitude 
of the A-doubling in this case is difficult due to uncertainty 
in the purity of the configurations. Qualitatively though, we 
would expect the effect of a 3 X~ state above the IIi state 
would be to push e-symmetry below f as observed. In ad- 
dition we observed A-doubling that was strongly vibrational- 
level dependent (3.69 x 10~ 4 cm -1 for the v = compared 
with 2.68 x 10~ 4 cm -1 for the v = 1), which is indicative of 
nearly-degenerate interacting states. We also measured a A- 
doubling term proportional to J 4 (denoted ko). This term can 
be explained by substituting AE — > AE + (Bn - B-£)J(J +1) 
in the A-doubling denominator; the first two terms in the Tay- 
lor expansion give the J(J +1) and the J 2 (J + l) 2 components. 
From this we can estimate that kp as 3k(Bn - B-£)/(En - E-£), 
which gives ko « 2 x 10~ 8 cm -1 . 



4. Theory 

In earlier theoretical study of HfF + fl9ll . the correlation cal- 
culations of the spectroscopic constants were performed in two 
different ways: in the first series 10 electrons from 5d, 6s shells 
of Hf and 2s, 2p of F were correlated while in the second series 
5s,5p outer core electrons of Hf and 1 s of F were also corre- 
lated. The rest of the 1 s - 4c/ inner-core electrons of Hf were 
excluded from the explicit treatment using the generalized rel- 
ativistic effective core potential (GRECP) method^ Q. While 
it was shown that the inclusion of the additional electrons in 
the 20-electrons case contributed significantly to the excitation 
energies, such 20-electron calculations were performed only 
for one internuclear Hf-F distance in a small space of many- 
electron basis functions due to computer limitations. This cal- 
culation was then applied for other distances as a "core correc- 
tion". In the present work a new 20-electron relativistic corre- 
lation calculation was carried out in order to consider electronic 
states that were not investigated in Jl9ll as well as to overcome 
the bottlenecks of the previous studies of HfF + . Eventually 23 
electronic states (with excitation energies up to 22000 cm -1 ) 
were considered in the present study instead of ten states in 



Some modifications were made in the computational pro- 
cedure to achieve better accuracy. We used a direct multi- 
reference configuration-interaction approach accounting for 
spin-orbit effects (SODCI) 1571 15811 as a method to treat 
both correlation and relativistic effects simultaneously. One- 
component basis functions (orbitals) are required in this method 
to construct many-electron spin-adapted functions (SAF's). In- 
stead of the orbitals obtained within the complete active space 
self-consistent field method (used in IH^P . the natural or- 
bitals of a one-electron density matrix averaged over the den- 
sity matrices of the states of interest were used. The latter 
density matrices were calculated within the scalar-relativistic 
coupled-clusters method with single and double cluster ampli- 
tudes (CCSD) using the cfour code f59ll . From numerical in- 
vestigation Il60ll it is known that the use of natural orbitals can 
provide faster convergence in terms of the number of SAF's re- 
quired to account for a given part of correlation energy. 

In the SODCI calculation the same atomic basis sets for Hf 



60-electron core GRECP for Hf is available at www.qchem.pnpi.spb.ru 



and F were used as in our previous study II 1911 . i.e., a generally- 
contracted basis set for Hf consisting of 6 s-, 5 p-, 5 d-, 3 f- 
and 1 g- type contracted functions (denoted as {6,5,5,3,1}) and 
an ANO-I {4,3,2,1} basis set for F 16 ill . However, a correction 
on the extension of the basis set was applied additionally in the 
present work. To evaluate the correction, two scalar-relativistic 
calculations were performed using the coupled clusters method 
with single, double and perturbative triple cluster amplitudes: 

(i) in the same basis set that was used at the SODCI stage, and 

(ii) with an extended basis set for Hf (produced by uncontract- 
ing d- and adding g-, h-, i- type basis functions) and an ANO-L 
{7,7,4,3} |6lll basis for F. 

To compute potential curves the described calculations were 
performed at 14 points in the range of 3.0 - 4.2 a.u. As shown in 
Table[3] the new calculations are in remarkable agreement with 
all of the experimental results. These improvements should 
enable more accurate calculations in other species with heavy 
atoms, where relativistic effects are extremely important. 

5. Outlook and Conclusions 

As shown in Figure[T]and discussed in the introduction, sev- 
eral different optical transitions may be necessary for the eEDM 
experiment. We have measured suitable candidates for all of 
these: (a) can be accomplished using the 3 no+ intermediate, (b) 
could use the 3 <E>2 level, and (c) could use one of the 3 no lev- 
els or the 3 <I>2 level. Another important value for the eEDM 
measurements is the size of the A-doubling in the 3 Ai J — 1 
level; it is the opposite parity levels in this state that are mixed 
in an electric field to polarize the molecule. The necessary elec- 
tric field for full polarization, G po i, is approximately where the 
Stark energy is larger than the energy difference between the 
two parity states: thus £ po / ss co e f/2nd m f J3l, where co e f is 
the A-doubling splitting (in angular frequency) and d m f is the 
molecular-frame electric dipole moment of the molecule (4.3 
Debye for HfF + ifiolo . We can use our measurement of 5a to 
estimate that a> e f — 2n x 45 a « 2n x 740 kHz, which is at 
least an order magnitude larger than previously predicted ifloTl . 
Nonetheless, this means that £> po i is still under 1 V/cm, which is 
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Figure 7: The eEDM induced Stark splitting ( AE) for / = 1 hyperfine levels of 
the 3 Aj state of HfF + : between F = 3/2 nip = ±3/2 levels (solid line), between 
F = 3/2 rtif = ±1/2 levels (doted line), and between F = 1/2 mp- = ±1/2 
(dashed line). 



spectrum II63I1 and a high repetition rate Ti:sapphire amplifier 
1 6411 . we hope to cover over 2000 crrT 1 with comb-vms. This 
ability to rapidly cover thousands of cirT 1 with high sensitivity 
and high resolution will be a powerful new tool for the study 
of ions and radicals for many applications. In addition, appli- 
cations to other many other spectral regions are possible by us- 
ing different comb sources combined with non-linear optics for 
covering the near-IR Jil S or mid-IR i3 H ^ III ■ 
These sources plus readout systems using different cameras or 
optical up-conversion 17211 create the possibility of ion spec- 
troscopy anywhere from the visible to the mid-IR. 
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important for the eEDM experiment as larger fields tend to lead 
to more issues with systematic errors. As can be seen in Fig- 
ure|7] numerical results for 3 Aj J — 1 hyperfine levels confirm 
that at 1 V/cm the splitting in the \nif \ =3/2 levels due to eEDM 
related Stark shift is almost saturated at 2d e & mt . This saturation 
occurs as the molecule becomes fully polarized. These results 
were obtained as described in [62]. 

We can also use the improved theory calculations to better 
estimate the radiative lifetime, another important consideration 
for the eEDM experiment. Using the measured energy separa- 
tion and calculated dipole moments, we estimate the lifetime of 
3 Ai v = to be about 2 s 11911 . This is encouraging, as it sets the 
ultimate limit on the coherence time achievable in the experi- 
ment. The high precision that we obtain for rotational constants 
should enable microwave spectroscopy with minimal search- 
ing, which will then be able to resolve the fine and hyperfine 
structure in the 3 Ai level. These measurements are important 
for estimating systematics and for checking the accuracy of the 
calculated effective electric field that the electron experiences 
(£;„,). We predict that the 'E + (v = 0) f = 1 <- J" = transi- 
tion in 180 HfF + will be at 18.290(2) GHz and the 3 Ai(v = 0) 
/' = 2 <— J" = 1 transition doublet will be at centered at 
35.869(2) GHz. To get these frequencies, we used the best mea- 
surements of B e and u) e from Table[3]to calculate D e = 4Bj/u%, 
then used the measured Bo and calculated D e to determine the 
microwave frequencies. 

We are currently working on improving the comb-vms sys- 
tem and applying it towards characterization of ThF + , which 
has several advantages over HfF + for the eEDM experiment. 
However for this species, the ab initio calculations are even 
more challenging, leading to larger theoretical uncertainties. 
The first spectroscopy of ThF + using PFI-ZEKE plus some LIF 
(unassigned) has recently been published 111 311 . but, as was the 
case with HfF + , the higher excited states remain uncharacter- 
ized. By using highly non-linear fiber to broaden the comb 
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Table 1: Fitted constants for observed transitions in 180 HfF + in cm 1 . Quoted errors are at 95% confidence level, statistical only. For bands recorded with 
the wavemeter uncalibrated we include an additional 0.006 cm -1 uncertainty in vo. Internal consistencies in the data suggest that for the Ai> = bands, which 
typically suffer from partially resolved isotope structure, systematic errors in the fit quantities may exceed the quoted statistical uncertainties in some cases by a 
factor of two or three. Values without uncertainties were fixed in the fits. All constants are given for the vibrational levels involved in the transition and are not 
equilibrium values. The A-doubling terms, k" , k' , and k' D are generic terms proportional to J(J + 1) in the ground and excited states and J 2 (J + l) 2 in the excited 
state respectively (see text for more details). 









B" 


B' 


D" 

[io- 7 ] 


D' 

[io- 7 ] 


k" 

[io- 4 ] 


k' 

[io- 4 ] 


[io- 9 ] 


'Ill* 


(0, 1) 


12217.369(2) 


0.30335(2) 


0.28115(3) 


1.88(8) 


1.81(8) 




3.69(2) 


9.7(7) 


'III* 


(1,2) 


12136.012(3) 


0.30180(5) 


0.27973(5) 


1.80(14) 


1.74(14) 




2.68(4) 


7.4(1.3) 




-'Z+ (0, 0) 


13002.189(12)+ 


0.30474(20) 


0.28104(20) 


1(2) 


1(2) 




3.55(4) 





3 IV 


-*£+ (2, 0) 


12304.400(3) 


0.30481(5) 


0.28096(5) 


1.78(12) 


1.78(12) 




-3.82(1) 





3 IV 


-'2+ (3, 1)° 


12216.901(10) 


0.30335* 


0.27958(20) 


1.8(3) 


1.8(3) 




-3.8(2) 





^ 


(0, 1) 


10109.877(6) 


0.30333(6) 


0.28382(6) 


1-9(5) 


1-9(5) 




-3.77(5) 





3 n 0+ 


^'S+ (0, 0) 


10401.723(13) + 


0.30500(10) 


0.28437(11) 


2.3(1.2) 


2.2(1.2) 








3 n 0+ 


(1,0) 


11114.653(7) 


0.30482(4) 


0.28281(4) 


1.71(14) 


1.72(13) 










^'1+ (0, 0) 


13254.302(7) 


0.30478(10) 


0.28967(9) 


1.7(5) 


1-9(5) 








3 5T ■ 


(1,0) 


13953.799(6) 


0.30483(5) 


0.28808(4) 


1.78(12) 


1.98(13) 








3 0> 2 < 


- 3 A, (0, 0) 


13933.340(12)+ 


0.29900(10) 


0.27720(10) 


2.2(3) 


2.1(3) 


t 


-0.01(3) 




3 o 2 < 


- 3 A! (1,0) 


14613.782(7) 


0.29888(5) 


0.27590(5) 


1.71(15) 


1.59(15) 


t 







3 o 2 < 


- 3 A! (2,1) 


14534.839(7) 


0.29745(8) 


0.27488(8) 


1.7(4) 


1.6(4) 


t 







3 n 0+ 


^ 3 Aj (1,0) 


10137.723(7) 


0.29887(5) 


0.28277(5) 


1.70(20) 


1.66(20) 


-0.14(2) 






3 n - 


^ 3 Aj (1,0) 


9948.624(6) 


0.29891(2) 


0.28272(2) 


1.85(5) 


1.81(4) 


-0.123(6) 








- 3 A 2 (0, 0) 


10852.757(13)+ 


0.29915(9) 


0.28102(9) 


1.7(4) 


1.7(4) 





3.59(3)* 




? Av 


= 


13729.918(17)+ 


0.29860(60) 


0.28680(60) 


3(3) 


3(3) 





-0.2(3) 




? Av 


= 


13822.270(14)+ 


0.29980(10) 


0.28860(10) 


1.9(2) 


2.2(2) 










? Av 


= 


14636.162(13)+ 


0.29620(30) 


0.27600(30) 


1.7(1.5) 


1.8(1.5) 











o3 rii <-'X+(3, 1) values assigned using a manual fit to multiple isotopes due to challenges of determining the line centers of the weaker lines in the dense 
spectrum. 

* Value fixed to the fitted value from the 1 TIi < — '2 + (0, 1) transition. 

+ Individual isotopes were not fully resolved, thus the fit was done to averaged line positions. We report isotope corrected values with the error due to the 
uncertainty in the isotope shifts. 

t 3 <I>2<— 3 Aj transitions exhibited line-doubling at high J", but the doublets were not resolved well enough to accurately determine k". 

* Each transition was a doublet with a splitting given by V X J' (J' + 1) due to A-doubling in ' IT i . 



Table 2: Derived constants for observed states in 180 HfF + in cm Values for To are given to the v = levels, with the X'S + v = level set to 0. Values for 
T e are to the minimum of the potential curves again withX'S + set to 0. AG 1/2 is only given for states where we had a direct measurement. All equilibrium values 
assume a Morse potential unless otherwise noted. Quoted uncertainties are 95% and are statistical, i.e., they do not account for any deviation from the Morse 
potential. ST e = T ls0 - T™ is an electronic isotope shift relative to 'S + for each state as discussed in the text. 





T„ 


AGi /2 


T 






LO e X e 


a e [10- 3 ] 




'z+ 





784.820(12) 





0.30558(3) 


790.76(11) 


2.97(5) 


1.50(2) 





3 A, 


976.930(10) 




991.83(74) 


0.29963(5) 


760.9(1.5) 


2.78(21 ) 


1.45(8) 


-0.049(12) 


3 A 2 


2149.432(18) 














-0.055(12) 


1 n 1 


13002.189(12) 




13046.04(21) 


0.28186(5) 


702.9(4) 


2.70(15) 


1.42(6) 


-0.071(6) 


3 n! 


10894.697(14) 




10933.77(6)+ 


0.28454(8) 


712.382(21)+ 


2.51(6)+ 


1.43(4) 


-0.054(7) 


3 n - 






10248.34(36)* 


0.28517(17)* 








-0.033(16) 


3 n 0+ 


10401.723(13) 


712.930(15) 


10437.44(36) 


0.28518(13) 


719.5(7) 


3.3(3) 


1.63(11) 


-0.047(7) 




13254.302(7) 


699.497(9) 


13296.93(28) 


0.29046(14) 


705.5(5) 


3.0(3) 


1.59(10) 


-0.106(7) 


3 <t>2 


14910.270(16) 


680.442(14) 


14963.15(33) 


0.27778(14) 


684.5(7) 


2.05(31) 


1.16(14) 


-0.044(16) 



+ Calculated without assuming a Morse potential. 

* Assuming oj e , aj e x e , and a e are the same as for the 3 fIo + . 
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Table 3: Summary of states. Comparison of derived molecular constants from this work and 
from the experiments of Barker et al. |[l2ll with the old theoretical calculations of Petrov et al. fl^l 
and the improved theory discussed here. The theoretical values of B e from tl9ll were computed 



from the equilibrium bond length. Experimental uncertainties are 


converted to 95% (2<r). 


State 


Constant 


This Work 


[12] 


U9J 


New theory 


ly+ 


T 

if 


U 


u 


u 


u 




B, 


0.30558(3) 


0.304(10) 


0.3082 


0.309 






790.76(11) 


791.2(1.0) 


751 


792 


3 A 

Al 


if 




yyj{Z ) 




izzy 




B c 


0.29963(5) 


0.301(10) 


0.2994 


0.301 




oj c 


760.9(1.5) 


761.3(2.0) 


718 


754 


3 A„ 


T 
if 


z m-y .t+jzy ioj 


9 1 c 1 7/"7nY'" 

zi j 1 . / yzxj) 


ZoU / 


HQA 




B, 


0.29915(9)* 


0.300(10) 


0.2997 


0.302 




co e 


- 


762.3(2.0) 


719 


766 


3 A 

A3 


if 




jyjlyZ ) 


A 11 A 


jyyj 




B, 




0.308(10) 


0.3004 


0.301 






- 


761.5(2.0) + 


721 


757 




T 
if 






1 1 J IV 






B e 


_ 


— 


0.2981 


0.298 




CO,. 


- 


- 


696 


747 


3 n„ 
iio- 


T 
if 


1 mzts ^zt("^^~\ 

f UZ'+oO^f^O ) 




1 1 qi n 
1 1 V1U 


1 DADO 




B, 


0.28517(17) 


_ 


0.2848 


0.286 




co c 


- 


- 


689 


716 


110+ 


T 






1Z 1 VD 


1UDJO 




Bf 


0.28518(13) 


_ 


0.2854 


0.285 




co e 


719.5(7) 


- 


699 


724 




T 
if 






IZOoO 


1 1 fKR 
1 lUJo 




Bf 


0.28454(8) 




0.2835 


0.285 




co c 


712.382(17) 




687 


712 


3 n 2 


Tf 






14438 


13452 




Bf 






0.2848 


0.287 




to e 






703 


745 




Tf 


13046.3(3) 




14784 


13493 




Bf 


0.28186(5) 




0.2805 


0.283 






702.9(4) 




679 


699 


3 sr 


Tf 


13296.93(28) 






13773 


Bf 


0.29046(14) 






0.292 




0> e 


705.5(5) 






716 



+ AGi /2 . 
"Bo- 
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Table 3: Summary of states, (continued) 



State 


Constant 


This Work 


[12J 


[19J 


New theory 


% 


T e 








14757 


B, 








0.292 




w e 


_ 


_ 


_ 


711 


J <1>2 


T e 


14963.15(33) 






15284 




B c 


0.27778(14) 






0.278 




k> e 


684.5(7) 


_ 




671 


G>3 


T, 








17457 




B, 








0.277 




(O e 


_ 


_ 




658 


r 4 


T e 








18312 




B e 








0.289 




w e 


_ 


_ 




641 




T e 








19167 




B e 








0.276 




w ( . 


_ 


_ 


_ 


691 


ni 


T, 








19332 




B, 








0.279 




<jj c 


_ 








698 




T e 








20074 




B ; , 








0.280 




w e 








748 




T, 








20330 




B, 








0.288 




w ( , 


_ 


_ 




610 


n 2 


T, 








20338 




B, 








0.277 






_ 


_ 


_ 


665 


3 4> 4 


T e 








20769 




B e 








0.283 




<y e 








740 




T e 








21415 


B, 








0.276 




w ( , 








665 


3y+ 


T e 


21200° 






21694 


B e 








0.277 




0» e 








658 



° Estimated from A-doubling in the 3 ITo (see text). 
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